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ABSTRACT 

It is found that there exists an empirical linear relation between 
the high frequency u higll and low frequency u\ ow of quasi-periodic os- 
cillations (QPOs) for black hole candidate (BHC), neutron star (NS) 
and white dwarf (WD) in the binary systems, which spans five or- 
ders of magnitude in frequency. For the NS Z (Atoll) sources, v high 
and v iow are identified as the lower kHz QPO frequency and horizon- 
tal branch oscillations (HBOs) z/hbo (broad noise components); for the 
black hole candidates and low-luminosity neutron stars, they are the 
QPOs and broad noise components at frequencies between 1 and 10 
Hz; for WDs, they are the "dwarf nova oscillations" (DNOs) and QPOs 
of cataclysmic variables (CVs). To interpret this relation, our model 
ascribes v high to the Alfven wave oscillation frequency at a preferred 
radius and v iow to the same mechanism at another radius. Then, we can 
obtain u\ ow = 0.08 u high and the relation between the upper kHz QPO 
frequency z/ 2 and HBO to be z/hbo — 56 (Hz) (z/ 2 /kHz) 2 , which are in 
accordance with the observed empirical relations. Furthermore, some 
implications of model are discussed, including why QPO frequencies 
of white dwarfs and neutron stars span five orders of magnitude in 
frequency. 

Subject headings: accretion — stars: neutron — white dwarfs — stars: oscillations - 
X-ray: binaries 
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1. Introduction 

To search the timing properties of the aperiodic variability of accreting X-ray binaries, 
the Rossi X-ray Timing Explorer (RXTE) was launched in the end of 1995, and has 
risen many discoveries (e.g., van der Klis 2000, 2005, 2006 for recent reviews). The quasi 
periodic oscillations (QPOs) and broad noise components are observed in the power spectra 
of neutron stars (NSs) and black hole candidates (BHCs) in low-mass X-ray binaries 
(LMXBs), which covers many QPO frequency bands, such as the kilohertz (kHz) QPO 
~ 200 - 1300 Hz that has been found in about 28 LMXBs (e.g., van der Klis 2006) and 
typically occurs in pairs: the upper kHz QPO u 2 and lower kHz QPO vi, the ~ 15 — 70 
Hz horizontal branch oscillations (HBOs) z^hbo i n Z sources or broad noise components in 
Atoll sources (on the definition of Atoll and Z, see Hasinger & van der Klis 1989; Hasinger 
1990). The numerous QPOs and broad noise components in the 0.1 - 1300 Hz range 
are found to follow the tight correlations covering three orders of magnitude in frequency 
(see, e.g., Psaltis, Belloni & van der Klis 1999, hereafter PBK99; Belloni, Psaltis & van der 
Klis 2002, hereafter BPK02). For the bright Z sources, the two QPO features are the HBO 
and the lower kHz QPO, respectively. For the BHCs and low-luminosity sources, they are 
the low-frequency QPO and a broad noise component at frequencies between 1 and 10 Hz. 
PBK99 and BPK02 have demonstrated that these NS and BHC low and high frequencies 
follow a tight correlation over 3 orders of magnitude in frequency; the correlation can be 
approximately described by a linear relation as u\ ow = 0.08u high (Titarchuk & Wood 
2002, hereafter TW02, see also figure 1). 

Moreover, dwarf nova oscillations (DNOs) and quasi-periodic oscillations (QPOs) were 
discovered a few decades ago (Warner & Robinson 1972; Patterson et al. 1977; Warner 
2004 for a review). Motivated by Warner & Woudt (2002ab), Mauche (2002) studies 
the data of the Extreme Ultraviolet Explorer deep survey photometer and the Chandra 
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X-ray Observatory Low Energy Transmission Grating Spectrograph to investigate the 
extreme-ultraviolet (EUV) and soft X-ray oscillations of the dwarf nova SS Cyg in outburst. 
To combine with the optical data of Woudt & Warner (2002) for VW Hyi, he extends the 
zAow-^high correlation (u\ ow « 0.08 fhigh) for NS/BHC LMXBs nearly two orders of magnitude 
in frequency. This correlation identifies the high-frequency quasi-coherent oscillations 
(so-called "dwarf nova oscillations" , DNOs) of cataclysmic variables (CVs) with the lower 
kHz QPOs of LMXBs, and the low-frequency QPOs of CVs with the HBOs of LMXBs. 

Furthermore, systematically studied and improved by Warner et al (2003) and Warner 
& Woudt (2004ab), the 27 white dwarf CVs and LMXBs are found to share a common 
correlation in the QPO frequencies over nearly five orders of magnitude in frequency, which 
also establishes the connection between the CV and LMXB QPOs. In addition to the 
frequencies, the DNOs of CVs and lower kHz QPOs of NS LMXBs are similar in that they 
have relatively high coherence and high amplitudes, and sometimes occur in pairs (Warner 
et al 2003; Warner & Woudt 2004ab; 2006; Pretorius et al. 2006). However, while the 
similar QPO phenomena occur in CVs and BH LMXBs, they are not totally identical (e.g. 
Kluzniak et al. 2005). 

On the origins for QPO productions, some models exist for LMXBs. The beat 
frequency models (e.g., Miller, Lamb & Psaltis 1998; Strohmayer et al. 1996; Alpar & 
Shaham 1985) take the NS spin-frequency as one of the characteristic frequencies, however 
it is commented recently that this model still confronts the severe arguments 
or modifications in interpreting the varied twin kHz QPO separations and the 
correlation between the upper and lower kHz QPO frequencies (see, e.g., Swank 
2004; Belloni et al. 2005, 2007; van der Klis 2006; Zhang et al. 2006). In the 
relativistic precession models (see, e.g., Stella & Vietri 1999; Stella et al. 1999; Morsink 
2000), all QPO frequencies arise from general relativistic frequencies in the accretion disk, 
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i.e., the upper and lower kHz QPOs (z/ 2 and u\) are ascribed to the orbital Keplerian 
frequency v K and the periastron precession of it at the same orbit, respectively, whereas the 
HBO frequency is ascribed to the nodal frequency of the tilt disk by the gravitomagnetism 
precession, frame-dragging effect of Einstein's general relativity or the Lense-Thirring effect 
(Stella & Vietri 1998; Stella et al. 1999). Although the empirical quadratic relation between 
the HBO and upper kHz QPO can be derived in the relativistic precession model, it needs 
the high spin frequency, such as ~ 900 Hz in central object, and large moment of inertia 
(see, e.g., Psaltis et al 1999), which are inadequate for WD sources because of their weak 
gravitational field and slow rotation (see, e.g., Mauche 2002; TW02; Warner & Woudt 
2004ab). 

More promising are models that interpret QPOs as manifestations of disk accretion 
onto any low-magnetic field compact objects, such as the transition layer model of Titarchuk 
and collaborators (Titarchuk, Lapidus & Muslimov 1998; Titarchuk et al 1999; Titarchuk 
& Osherovich 1999, 2001; Titarchuk et al 2001; TW02). In the model of TW02 the high 
frequency z/Mgh (identified as the lower kHz QPO frequency v\) is ascribed to the Keplerian 
frequency u K at the outer (adjustment) radius and the low frequency z/i ow (identified as 
the HBO) is ascribed to the magnetoacoustic (MA) oscillation frequency z/ M a, and they 
demonstrate that the observed correlation of the low-high frequency is perfectly described 
by the dependence of the inferred MA frequency %a on the Keplerian frequency v K . 
Moreover, as for the QPO models that are nothing to do with the presence or absence 
of the hard surface of the compact object, there exist the mode oscillation models in the 
accretion disc (see, e.g., Li & Narayan 2004; Rezzolla, et al 2003ab; Kato 2001; Nowak et 
al 1997; Wagoner 2001, references therein) and the orbital resonance model (Abramowicz 
et al. 2001, 2002, 2003; Rebusco 2004; Petri 2005; Horak k Karas 2006), which needs the 
further refinements if applied to interpret the low and high frequency correlations. 
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In this paper, based on the Alfven wave oscillation model that has successfully explained 
the kHz QPOs of NS LMXBs (Zhang 2004), we seek to offer a unified prescription to explain 
the similar correlation between the low and high frequencies found in CVs and LMXBs. 
In §2, our model formulation is described, and the derivation of the low-high frequency 
correlation is presented as well. The consequences and conclusions are summarized in §3. 

2. Description of the model and its applications 

In the model by Zhang (2004), it ascribes the upper and lower kHz QPO frequencies 
[y 2 and ui) to the Keplerian frequency and the MHD Alfven wave oscillation frequency at a 
special radius r, where 

v 2 {X) = 1850AX 3/2 (Hz) = 1295(A/0.7)X 3/2 (Hz) , (1) 

with X=R/r, A = (m/R^) 1 / 2 , the scaled stellar radius R 6 = R/10 6 cm and mass in the unit 
of solar mass m = M/M , 

Ul (X) = u 2 (X)X(l - Vl-X) 1/2 = 1295(A/0.7)F(X) , (2) 

F(X) =x 3 /(i + VT zr x) 1/2 ■ (3) 

The inferred averaged value of the parameter X (X) ~ 0.88 comes from the best 
fit of the kHz QPOs of four detected typical kHz QPO sources and it is not 
an optimization, then from (X) ~ 0.88 it infers that the position of producing 
the kHz QPO is close to NS surface, about ~ 10% of star radius away (Zhang 
2004). It is remarked that our model (Zhang 2004), successfully applied in interpreting 
the kHz QPO correlations of NS LMXBs, just needs the accretion flow in the disk around a 
gravitational source with the dipole magnetic moment inside an innermost boundary, and 
it does not necessarily require the presence or absence of a hard surface of compact object. 



- 7- 



To study the mechanisms of the low and high frequencies detected in 
NS/BH/WD systems by means of our kHz QPO model (Zhang 2004), we 
ascribe the high frequency i/^igh to the Alfven wave oscillation frequency at 
a preferred radius r and the low frequency u\ ow to be the same mechanism 
occurred at another radius of about <pr, where is a parameter determined by 
the observation. Therefore, these low and high frequencies (i^ ow and z/high) are 
conveniently arranged as follows 

^high = Vi{X) , (4) 

^iow = vi{(j)~ l X) , (5) 
and then we have the ratio between the low and high frequencies, 

2Z-=W,X), (6) 

^high 

where ratio quantity K(<f>,X) is defined as 

K( ( p,X) = 0.0M(^)- 3 g( ( p,X), (7) 



9^X) = [-1±^^=]^<1. (8) 
1 + y/l-^X 

Because the position parameter X=R/r is in the range of < X < 1 from its 
definition, we find that K((f>, X) is mainly determined by <fi but it is insensitive 
to X, and even more important the ratio quantity K(<p, X) is independent of the 
parameter A, where A 2 = m/i?| represents a measurement of the average mass 
density of compact object. Therefore Eq.Q can be applicable in any compact objects, 
from WDs to NSs. For the possible values of = 2.2 and X=0.88, from Eq.© and Eq.® 
we obtain the ratio K ~ 0.08, which satisfies the need of empirical relation by TW02. 

The correlation f low versus i^igh is plotted in figure 1, together with the measured data 
of CVs and LMXB sources from BPK02, Mauche (2002) and Warner & Woudt (2004), 
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respectively. The three fitted lines, from top to bottom, represent v\ow/ ^high = 0.11, 0.08, 0.06, 
which can be obtained by setting the parameter values of <fi = 2.1,2.2,2.3 and X=0.88, 
respectively. It is noticed that the agreement between the model and the observed QPO 
data is quite good for the convenient choice of parameters X and <fi. It is stressed that the 
ratio K((f>,X) in Eq.([7|) has nothing to do with the parameter A and is almost independent 
of the parameter X. Therefore, the z/ low versus z/ high plot in figure 1 will represent the 
properties of any compact objects (CVs and LMXBs) that reflect a common feature of the 
accretion flow around the gravitational sources and have no direct correlations to the mass, 
radius, hard surface, magnetic field and spin frequency of the specific source. 

If we assume the position parameter X to be similar for the QPO sources of WDs 
and NSs, then for the typical choice of mass and radius of WD and NS (e.g., Shapiro 
& Teukolsky 1983), such as M WD = 1.0M o , R WD = 5 x 10 8 cm, M NS = 1.4M and 
R NS = 1.5 x 10 6 cm, respectively, we can obtain, 

^high(WD) _ A WD F(X WD ) _ M WD R WD 3 / 2 F(X WD ) 4 

^igh(NS) ~ A NS F(X NS ) ~ 1 M NS ' 1 R NS ' F(X NS ) ~ U 

where we have set the position parameter Xwd ~ Xns or wj^j ~ 1 (see Eq. EJ. If the 
WD mass (radius) is lower (higher) than the typical value or the condition Xwd > X^s 
happens, then the above ratio in Eq.Q can be even as low as ~ 10~ 5 , which may interpret 
the extended low-high frequency correlations in CVs and LMXB systems in five orders 
of magnitude in frequency, as shown in figure 1. In fact, there are two parameters to 
modulate the QPO frequency in Eq.(j2j) or Eq.(j4j), A and X, representing the compactness 
of source and the position to exhibit QPO respectively. Therefore for the same source (A is 
same), the parameter X will modulate the QPO frequency, while for the different sources 
the parameter A will modulate the QPO frequency. However, the fact that the ratio of 
low frequency to high frequency is independent of parameter A and weakly depends on 
the parameter X can interpret the QPO similarities in different sources and in different 
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accretion rates of the same source (see figure 1). 

As an assumption that u\ ow is identified as the HBO frequency ^hbo f° r Z sources (or 
Broad band noises for Atoll sources), i.e., i^hbo = iAaw, we obtain the theoretical correlation 
between the HBO frequency ^hbo an d the upper kHz QPO frequency vi from Eqs.fll]) and 
(JSD as, 

z/hbo = 74 (Hz) (^)- 3 (A/O.7)- 1 (z/ 2 /kHz) 2 /(0A) , (10) 

/(0A) = 1/(1 + V ^A70) 1 / 2 . (11) 

After substituting = 2.2, the averaged values (X)=0.88 and A = 0.7 inferred from the 
typical kHz QPO sources (Zhang 2004), into Eq. fTOjl . we obtain the simplified theoretical 
correlation between i^hbo and v-i as 

z/hbo ^ 56 (Hz) (z/ 2 /kHz) 2 , (12) 

which is consistent with the observed empirical relation (see, e.g., Psaltis et al 1999; Swank 
2004; van der Klis 2000, 2006). 

3. Consequences and Discussions 

In this paper, a unified interpretation of QPO low-high frequency correlations observed 
in CVs and LMXB sources has been proposed by means of the Alfven wave oscillation 
model for kHz QPOs (Zhang 2004), where we identify the high frequency to be the 
lower kHz QPO frequency v\ and the low frequency ui ow to be the same mechanism of V\ 
at the different radius. While, the linear relation between v\ ow and fhigh can be obtained as 
^iow = 0.08z/hi g h, which is consistent with the detected empirical relation (see, e.g., TW02; 
van der Klis 2000, 2006; Psaltis et al 1999; PBK99; BPK02), firstly established by PBK99 
and BPK02 for LMXBs and subsequently extended by Warner & Woudt (2002ab), Mauche 
(2002), Warner et al (2003) and Warner & Woudt (2004ab) for WD binaries. In addition, 
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the quadratic correlation between the low frequency z/[ ow (or HBO ^hbo) and upper kHz 
QPO frequency u 2 is approximately obtained as ^hbo^Iow) — 56 (Hz) (i/ 2 /kHz) 2 , which is 
also in accordance with the observed empirical relation in LMXBs (see, e.g, van der Klis 
2000, 2006; Psaltis et al 1999; Stella et al 1999). 

Undoubtedly, these studies suggest that the same types of variabilities occurred in 
CVs and LMXBs can severely constrain the theoretical models (see, e.g., TW02; Mauche 
2002) and provide clues to understanding the nature of QPO phenomena. Apparently, a 
common QPO feature for a wide class of CV and LMXB systems seems to hint that the 
mechanism to produce the low-high frequency correlation has no strong direct dependence 
on the specific parameters of sources, such as the mass, radius, spin, presence or absence of 
a hard surface of compact object and magnetic field (see, e.g., Mauche 2002; TW02; Warner 
& Woudt 2004ab). 

In conclusion, the extension in frequency for five orders of magnitude in z/iow-^high 
relation among CV and LMXB systems may imply that the frequencies z/i ow and ^hi g h are 
intimately related to the parameter A = (m/i?g) 1//2 or the Keplerian frequency uk where 
the parameter A appears (see Eq.flT])) because the typical values of masses and radii of WDs 
and NSs can arise the variation in the parameter A for five orders of magnitude. However, 
clearly seen in our model, ui ow and z/hi g h satisfy the above condition and the ratio between 
them is almost a constant that is nothing to do with the parameter A. Furthermore, for the 
different CV and LMXB sources the variation of z/i ow or t^igh is modulated by the variation 
of parameter A, and then for the same source it is modulated by the position parameter X 
which depends on the accretion rate. This fact can explain why the QPO data distribution 
of CV and LMXB sources in tVv-^high plot is homogeneous for the different sources and 
for the different accretion rates of the same source (see figure 1). By the way, the QPO 
similarity phenomena are not special for the CV and LMXB sources, which have also been 
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found in micro- quasars and quasars with 3:2 ratio in their twin-peaked QPO frequencies 
(see, e.g., Abramowicz et al 2004), and they may reflect the fundamental characteristics of 
accretion flow around any astrophysical objects in universe. 

However, the kw -i/ high correlation of black-hole timing features are in 
between those of WD and NS (See Fig. 1). To interpret this, we mention 
the fact that for black-hole systems (as well as the low-luminosity neutron 
star-systems), u higll corresponds to a very broad feature rather than a narrow 
QPO (BPK02). Moreover, to successfully apply our model to BHC, which has 
a different environment from those of WD and NS, such as the hard surface and 
magnetic field in WD or NS, we are inclined to assume that a dipole magnetic 
field is rooted in the accretion disk of BH and origins at its innermost stable 
circular orbit (ISCO). So, by conveniently choices of the parameters like what 
have been done for WD and NS, the special radii can be constructed to account 
for the detected relation between z/ low and z/ high . 

Finally, we stress that the low and high QPO behaviors in CV and LMXB sources are 
so complex that no simple explanations have been likely to be adequate (see, e.g., van der 
Klis 2006). As declaimed by Warner and Woudt (2004ab), the existence of magnetically 
controlled accretion is in principle testable by the further QPO detections. Upon the QPO 
correlation in CV and LMXB sources, it essentially implies that the QPO phenomena 
generate in the accretion disc around any low-magnetic field compact objects, and therefore 
our model would provide a competitive choice. 

We thank T. Belloni, M. Mendez, D. Psaltis, C.W. Mauche, B. Warner and PA. Woudt 
for kindly providing the QPO data. Discussions with J. Petri, P. Rebusco, J. Hordk, V. 
Karas, G. Hasinger, T.P. Li, J. L. Qu, L.M. Song and S.N. Zhang are highly appreciated. 
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Fig. 1. — The low QPO frequency versus high frequency diagram for the LMXBs (filled 
squares: black hole binaries, open squares: neutron star binaries), 27 CVs (filled circles) 
plotted only once for each CV, and SS Cyg and VW Hyi (open circles). The QPO data of 
LMXBs are from Belloni, Psaltis & van der Klis (2002) and kindly provided by T. Belloni. 
The CV data are from Warner & Woudt (2004ab) and Mauche (2002), which are kindly 
provided by B. Warner & P. Woudt and C. Mauche. The three straight lines represent the 
theoretical curves fiow/^high = 0.11, 0.08, 0.06, respectively, from top to bottom. 



